Epstein-Barr Virus DNAemia and post-transplant lymphoproliferative disorder in pediatric solid organ transplant recipients

Background Pediatric solid organ transplant (SOT) recipients commonly have Epstein-Barr virus (EBV) DNAemia and are at risk of developing post-transplant lymphoproliferative disorder (PTLD). EBV DNAemia has not been analyzed on a continuous scale in this population. Methods All children ≤ 18 years of age who underwent SOT at a single center between January 1, 2007 and July 31, 2018 were included in this retrospective study. Transplant episodes in which PTLD occurred were compared to transplant episodes without PTLD. Multivariable logistic regression was used to identify factors associated with the development of EBV DNAemia and maximum height of EBV DNAemia. A Cox proportional hazards model was used to calculate hazard ratios for time to PTLD. Results Of 275 total transplant recipients and 294 transplant episodes, there were 14 episodes of PTLD. Intestinal and multivisceral transplant were strongly associated with PTLD (p = 0.002). Risk factors for the development of EBV DNAemia include donor and recipient positive EBV serologies (p = 0.001) and older age (p = 0.001). Maximum level of EBV DNAemia was significantly associated with development of PTLD (p<0.0001). Every one log (log10) increase in the maximum level of EBV DNAemia was associated with a more than doubling of the hazard on developing PTLD (HR: 2.18, 95% CI 1.19–3.99). Conclusions Transplant type was strongly associated with development of PTLD in pediatric SOT recipients. EBV serologies and age were associated with the development of EBV DNAemia and height of DNAemia. High levels of EBV DNAemia were strongly associated with an increased hazard for PTLD.


Introduction
Post-transplant lymphoproliferative disorder (PTLD) is a significant cause of morbidity and mortality in pediatric solid organ transplantation (SOT) recipients. Epstein-Barr Virus (EBV) DNAemia plays a key role in the pathophysiology of EBV-related PTLD but characteristics of pediatric SOT recipients who develop EBV DNAemia and patterns of EBV DNAemia have not been well defined [1]. Many centers monitor EBV DNAemia with established protocols for interventional measures such as decreasing immunosuppression [2]. These protocols may differ between centers and between different types of organ transplant recipients [3]. More data is needed to guide clinical decisions in specific scenarios.
There have been many studies on risk factors for PTLD in SOT recipients [1][2][3][4]. Overall, risk factors for development of PTLD within the first 12 months of transplant (early PTLD) include: type of organ transplant, T-cell depleting therapies such as anti-thymocyte globulin (ATG), and young age. Early PTLD is also strongly influenced by donor and recipient serologies as donor-derived EBV infection is a strong risk factor for PTLD. Risk factors for development of PTLD after 12 months post-transplant (late PTLD) include length and duration of immunosuppression, and type of organ transplant. Intestinal transplant recipients have been found to have the highest risk for PTLD, followed by lung, then heart, then liver, and lastly, kidney transplant recipients.
There is a higher incidence of PTLD in pediatric SOT recipients when compared to adult SOT recipients [1]. One reason is that children are at higher risk for primary infection, and primary EBV infection has a higher risk of causing PTLD [4]. One pediatric SOT study found that EBV viral load was significantly higher in cases of PTLD, when matched to controls [5]. Another study reviewing only pediatric patients found the following risk factors to be significant in a univariate model: age at transplant, use of basiliximab, steroid intensification, ATG, median peak EBV value level, and chronically elevated EBV DNAemia [6]. However, in the multivariable model, only age, use of steroids, and peak EBV value were significant.
Many previous studies in pediatric solid organ transplants have categorized EBV DNAemia into three categories: 1) chronic high load (CHL) carriers, 2) chronic low viral load carriers, and 3) those with sustained undetectable EBV levels, with chronic defined as a time period > 6 months [7][8][9][10]. In pediatric heart transplant recipients, CHL carriers are predisposed to PTLD, and pediatric intestinal transplant recipients are at intermediate risk [7,8]. In the pediatric liver and kidney populations, CHL carriers are not at higher risk of PTLD [9,10]. Despite low sample sizes within studies, this difference suggests that transplant type and immunosuppression play a large role in the development of PTLD.
There have been many studies on trying to categorize cutoffs for EBV DNAemia and risk for PTLD, as well as establishing a predictive model for PTLD [6,11,12]. While many of these studies have tried to categorize a certain EBV quantitative cutoff level for when patients are at risk for PTLD, there are confounding risk factors such as patient age, type of transplant and type of immunosuppression and not all models have been well validated. In addition, there is controversy between the use of whole blood vs plasma samples, and there is interlaboratory variability for quantification of EBV levels from the same samples [13].
The epidemiology of PTLD has been changing in recent years [14,15]. In addition, EBV monitoring has provided a wealth of data that has not been studied in a detailed manner.
While EBV DNAemia has been categorized in previous studies, it has not yet been studied on a continuous scale and available pediatric literature has not linked maximum height of EBV DNAemia to PTLD on a granular level. Therefore, we undertook a comprehensive retrospective descriptive study at a large pediatric SOT center to specifically evaluate risk factors associated with EBV DNAemia. We also wanted to investigate the predictive value of continuous EBV monitoring for PTLD in order to demonstrate its ongoing utility.

Material and methods
We performed a retrospective study of patients � 18 years of age who received a solid organ transplant at Duke University Medical Center between January 1, 2007 and July 31, 2018. Events were recorded by transplant episode. The population included heart, kidney, liver, lung, intestinal, and multivisceral transplant recipients. Characterization of the pediatric SOT population was performed using the United Network for Organ Sharing (UNOS) database, and data was also extracted from the center's electronic health records using the program Duke Enterprise Data Unified Content Explorer (DEDUCE) [16]. The diagnosis of PTLD and lymphoma were identified by diagnosis code. The diagnoses were next verified through chart review of pathology results and progress notes in the electronic health record. The study was approved and granted a waiver of informed consent by the Duke University institutional review board (IRB). Patient health information and identifiers were available within a secure workspace, but only de-identified, anonymized data can be exported out of this workspace.
At our center, each pediatric solid organ transplant program has developed their own methods of screening for EBV DNAemia. However, in general, quantitative whole blood EBV polymerase chain reaction (PCR) studies are done monthly for the first year post-transplant, then spaced out to every 3 months for the second and third year post-transplant, then either every 6 months or yearly afterwards if there are no issues. If EBV DNAemia was detected, quantitative PCRs were obtained more frequently. All EBV PCRs included in this study were done at the central laboratory of our center, which used the same machine over the study period, the ABI PRISM 7500/7500DX Sequence Detection System (Applied Biosystems) and Qiagen EBV PCR (ASR) reagents.
EBV DNAemia values were extracted from the electronic health record, including accompanying chart review for missing values. For most of the study period, the quantitative unit of EBV DNAemia was reported as copies/μL. Starting in mid-2017, this was changed to international units/mL (IU/mL), and a conversion factor of 1 copy/μL = 113.6 IU/mL was used to normalize these values in concordance with the previous quantification. The laboratory checks this conversion every 6 months and the conversion factor has not changed over the duration of the study.
Patients who developed EBV DNAemia after transplant were compared to patients who did not develop DNAemia. EBV DNAemia was defined as any detection of EBV virus by quantitative polymerase reaction (PCR). EBV donor and recipient serologies were based on what was entered into the patient's chart at time of transplant. Demographic and clinical characteristics were compared between the two groups, using Chi square and Fisher's exact tests for categorical factors and Mann-Whitney U tests for continuous factors. A logistic regression was done to calculate odds of EBV DNAemia using the variables of age, race, type of transplant, type of induction immunosuppression, and donor and recipient EBV serologies. These variables were either significant in the unadjusted model or risk factors that were denoted to be clinically significant a priori based on current knowledge.
The diagnosis of PTLD was confirmed by chart review. All cases of PTLD had accompanying histologic diagnoses which included EBER staining and description of atypical cells seen on pathology. Sources of biopsy tissue included lymph nodes, tonsils and adenoids, and the transplanted organ.
For the analysis of the height of EBV DNAemia, the maximum single EBV value was identified. That value was then log 10 transformed and an analysis of covariance (ANCOVA) was used to evaluate associations between the maximum EBV value and demographic and clinical factors. After categorization and taking the average of the log 10 transformed values, the data was back transformed, establishing the geometric mean. Logarithmic comparison of quantitative EBV values has been well described in the literature [17,18].
A Cox proportional hazards model was used to link clinical risk factors to time to PTLD. The event of PTLD, death, and the event of a subsequent transplant were used as censoring criteria. The last observed value was carried forward. A robust variance (sandwich) estimator was used to control for the fact that some children were included more than once, due to multiple transplants. Adjusted models incorporated variables that were significant in the unadjusted model as well as risk factors that were denoted to be clinically significant a priori.
All children who developed PTLD were EBV positive by quantitative PCR, and thus, only EBV positive children were included in following analyses ( Table 2). Fourteen children developed PTLD after transplant, for an overall incidence rate of 4.8%. Transplant type was significantly associated with PTLD (p = 0.001). Intestinal and multivisceral transplant recipients accounted for 21% and 14% of PTLD cases respectively, while only accounting for 2% and 3% of the total population that was EBV positive. There was also an association between PTLD and race (p = 0.04) as there were a total of six Asian children in the study who were EBV positive, and two developed PTLD. Asian children comprised 14% of the PTLD cases, but only 2.2% of the study population. Induction immunosuppression and age were not associated with PTLD in our study (p = 0.18 and p = 0.17, respectively). For the nine patients who developed PTLD within the first year, EBV values were plotted over time (Fig 1). For most patients, the general pattern was a steady increase in EBV values, followed by diagnosis of PTLD, and then a rapid decline.
Risk factors for EBV DNAemia were examined (Table 3). Significant risk factors for EBV DNAemia included: EBV serology at the time of transplant, age, and type of transplant in the unadjusted model, and only EBV serology and age in the adjusted model (Table 4). Heart transplant recipients were at lower odds for EBV DNAemia, with OR 0.44 (95% CI 0.22-0.89). Lung transplant recipients were similarly at lower odds, but sample size was limited. Donor positive/recipient negative (D+/R-) and donor positive/recipient positive (D+/R+) transplant recipients were at increased odds of developing EBV DNAemia compared to donor negative/ recipient negative (D-/R-) transplant recipients, with OR 3.90 (95% CI 1.55-9.80), and OR 4.83 (95% CI 2.02-11.55), respectively. Older children were also associated with an increased risk, with every one year increase in age associated with an OR of 1.10 (1.04-1. 16).
An analysis of the geometric mean of the highest value of EBV DNAemia ( A final analysis was performed evaluating the association between risk factors and time to PTLD, and the adjusted analysis is included (Table 6). Every log 10 increase in maximum EBV value was associated with a hazards ratio (HR) of 2.18 for development of PTLD (95% CI  (Tables 7 and 8). The trends in these tables show that multivisceral and intestinal transplant recipients were tested the most frequently. The number of transplant recipients who were EBV positive by PCR post-transplant did not seem to change over time. In addition, while the number of tests for liver transplant recipients during the first year decreased over the study time period, the number of tests done during the first year for non-liver transplant recipients (heart, kidney, lung, intestinal, and multivisceral transplants) combined did increase overall which suggests better adherence to EBV PCR monitoring over time.

Discussion
This is the first study to fully describe characteristics of post-transplant EBV DNAemia in the pediatric solid organ transplant population. In a raw analysis of 275 SOT pediatric recipients, type of transplant and race were significantly associated with higher risk for PTLD. Intestinal and multivisceral transplant recipients had the highest risk of developing PTLD, which is  supported by the literature. We interpret the findings of race with caution, as the overall sample size is small, and only two out of four Asians with EBV DNAemia developed PTLD. We were unable to find any other risk factors that were significantly associated with PTLD in our population, including type of induction immunosuppression, EBV donor and recipient serologies, or age. This may have been due to distinctions within our population or to limited sample size.
When we examine risk factors for EBV DNAemia, and for magnitude of EBV DNAemia, some common themes emerge. For EBV DNAemia, EBV serologies of the donor and the

PLOS ONE
recipient are an obvious risk factor. Age is also a practical risk factor as the risk of being exposed to EBV increases with age. Type of transplant is not associated with risk of EBV DNAemia in the adjusted model, but within the category, we see that EBV DNAemia tended to be more common in liver transplant recipients. Age was the most resilient predictor for the maximum height of EBV DNAemia, as it changed the least from the unadjusted to the adjusted model. Intestinal transplants did not have a significantly different height of EBV DNAemia, yet they still had a higher risk association with PTLD, which suggests that other factors such as length and duration of immunosuppression should be considered. Ultimately, those who developed PTLD had significantly higher maximum levels of EBV DNAemia, which has been seen in previous studies. 5 However, stratifying by type of organ transplant does matter; liver transplant recipients were more likely to have a higher maximum level of DNAemia compared to heart and intestinal transplant recipients but were not more likely to develop PTLD. This may be due to immunosuppression or organ-specific risk factors for PTLD. Interestingly, the association of EBV serostatus D+/Ris often associated with the highest risk for PTLD, and our findings in Tables 2 and 6 do not reinforce this. This may have been limited by sample size, or other unmeasurable factors in our population.
Our final analysis is significant in that higher levels of EBV DNAemia are associated with decreased time to PTLD, even when adjusted for other risk factors. It is also important to note that an earlier date of transplant was associated with a higher hazards ratio for PTLD. This may have been due to better EBV monitoring over time and the development of protocols which have standardized practice at our center. The significance of a higher hazards ratio for kidney transplants is not clear at this time. We did have much more liver transplant recipients and less kidney transplant recipients which is a limitation of this study and this may have affected the hazards ratio. Lung transplants did not have a high incidence of DNAemia as indicated by Table 3, but did have high hazard ratios for PTLD. This again may be due to a small sample size of lung transplant recipients within our population. In addition, many lung transplant recipients from our center are referred back to more local home institutions and EBV PCRs may be obtained externally. Those who stay at our center often have many complications such as rejection, and of course, the development of PTLD. Lastly, type of induction immunosuppression may be important given the trend towards decreased time to PTLD with the use of alemtuzumab, ATG, or high dose steroids. These findings may help with stratifying patients who may be at higher risk for PTLD. It is not immediately quite clear why risk factors associated with PTLD (Table 2) were different than risk factors associated with increased hazards of PTLD (Table 6). For the analysis of hazard ratios, we did add an extra variable of the maximum level of EBV DNAemia and the Cox regression analysis was adjusted for different variables, and this may have changed the analysis. There may also be statistical differences when examining development of PTLD vs time to PTLD.
Risk factors associated with increased odds of EBV DNAemia were slightly different than risk factors associated with increased hazards of PTLD. This may reflect the fact that while EBV DNAemia is a risk factor for PTLD, how clinical providers manage the EBV DNAemia may influence the development of PTLD. Also, for example, liver transplant recipients were more likely to have EBV DNAemia, but they are less likely to develop PTLD so other factors such as type of organ transplant and immunosuppression may be more important in the development of PTLD.
All pediatric SOT studies for PTLD and EBV are limited by small sample sizes and this single center study is no exception. The overall incidence of PTLD was relatively low. This may have be due to the quick responsiveness of clinicians to high EBV levels. One major limitation is that we were unable to adjust for time-varying exposures to immunosuppressive medications. Higher periods of immunosuppression could have led to increased incidence of PTLD, EBV DNAemia, and higher maximum EBV levels. On the other hand, increased EBV levels may have led clinical providers to decrease immunosuppression, which could have prevented PTLD. However, there is no clear way to delineate all the changes to maintenance immunosuppression as changes occur often and it is very difficult to quantitate such a time-varying risk factor in our analysis. We did have access to a wealth of data due to frequent EBV monitoring and were able to examine multiple risk factors for EBV DNAemia and the maximum height of EBV DNAemia. However, this study was not adequately powered to examine the numerous interactions between risk factors that may occur in predicting PTLD. Another issue may be the introduction of testing by indication bias, where transplant recipients who had positive and higher levels of EBV were tested more frequently than others. We did note that for the most part, while those who did test positive for EBV by PCR were more likely to receive a higher number of total tests, those who tested negative for EBV by PCR also received a fair number of tests, which denotes that teams were adhering to their protocols (S1 Data).
Other limitations include identifying PTLD and lymphoma cases by diagnosis code and while we did verify each diagnosis, we may have undercounted the number of cases as we may have missed some cases if they were not coded a certain way. Our center uses whole blood sampling for EBV which is much more sensitive but not as specific. While other centers may have moved towards plasma sampling, the more sensitive measure may have led to less cases of PTLD. We did not establish specific cutoffs for the risk of PTLD, but our findings have been shown in a way that is much more useful to individual centers as our final analysis looks at the risk of PTLD associated with each change in log of the quantitative level. We only included EBV quantitative PCRs done at the central laboratory at our center, and we may have missed values that were obtained externally. The decision was made to not include external values in this study as that may have introduced interlaboratory variability to this study. In addition, EBV results from these external laboratories would have been very hard to find in our electronic records, as our center utilized paper charting prior to 2014 and they may not have made their way into each scanned chart from that era.
Despite the study limitations, we did gain some valuable insights into the nature of EBV DNAemia and PTLD in our study population. Further prospective or multicenter studies are needed to link significant clinical factors and EBV DNAemia with PTLD and other clinical outcomes.

Conclusion
Type of transplant and age were associated with PTLD in pediatric SOT recipients. Age and EBV donor and recipient serology remain key risk factors in the development of EBV DNAemia and the maximum height of EBV DNAemia. Maximum level of EBV DNAemia, and year of transplant were associated with time to PTLD in our population, and there was a trend towards type of induction and time to PTLD. While this study has shed some light on risk factors for EBV DNAemia in pediatric SOT recipients, further studies are needed in order to fully characterize the relationship between EBV levels and PTLD.